Quantum superposition of three macroscopic states 
and superconducting qutrit detector 
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Superconducting quantum coherent circuits have opened up a novel area of fundamental low- 
temperature science since they could potentially be the element base for future quantum computers. 
Here we report a quasi-three- level coherent system, the so-called superconducting qutrit, which has 
some advantages over a two-level information cell (qubit), and is based on the qutrit readout circuit 
intended to measure individually the states of each qubit in a quantum computer. The designed and 
implemented radio-frequency superconducting qutrit detector (rf SQUTRID) with atomic-size ScS- 
type contact utilizes the coherent-state superposition in the three- well potential with energy splitting 
AEoi/kB ~ 1.5 K at the 30th quantized energy level with good isolation from the electromagnetic 
environment. The reason why large values of A_Eoi (and thus using atomic-size Nb-Nb contact) 
are required is to ensure an adiabatic limit for the quantum dynamics of magnetic flux in the rf 
SQUTRID. 

PACS numbers: 03.75.Lm, 74.50.-|-r, 85.25.Cp 



I. INTRODUCTION 

The phenomenon of the superposition of states of 
a macroscopic object predicted for a superconduct- 
ing quantum interferometer device (SQUID) in the 
low-dissipation limiti^— was revealed in spectroscopic 
experiments 1^1^ Although much progress has been made 
in demonstrating the coherent quantum behavior of var- 
ious systems with Josephson junctions,^ there has not 
been an experimental presentation of a readout device 
based on the quantum superposition of macroscopically 
distinct states in flux qubits. The dependence of this 
fundamental property of quantum mechanics, the super- 
position of states, on the symmetry of the potential en- 
ergy in flux qubits can be taken as a basis for creat- 
ing a radio-frequency superconducting qubit detector (rf 
SQUBID) similar to the manner of how the Josephson 
current-phase relation is used in building rf SQUIDs.'' 
This device would be a principal element in quantum 
readout circuits meant for weak continuous measurement 
of states of the flux qubits incorporated in the quantum 
computer architecture. 

However, there is a pitfall on this way. In the litera- 
ture, the flux qubit (a superconducting ring closed by 
a superconductor-insulator-superconductor (SIS) junc- 
tion) is described as the superposition of the two states 
l^i) = ci(t) + C2(t) 1*2) that appear in this quan- 
tum system with double-well potential. If the tunnel- 
ing amplitude is large (in other words, if coeSicients ci 
and C2 vary quickly enough), the system behavior be- 
comes adiabatic when changing the external magnetic 
flux; that is, it can be considered in terms of quasi- 
stationary superposition levels. In this case, it becomes 
impossible to distinguish between the experimental char- 
acteristics of a common classical SQUID in a nonhys- 
teretic regime and a double-well SQUBID. Both devices 



behave as parametric inductances (Josephson inductance 
for SQUID and quantum inductance for SQUBID), with 
both inductances being maximized at the same exter- 
nal flux, (f>e = (n + l/2)$o (where n is an integer and 
$0 = h/2e is the superconducting flux quantum), so that 
some additional evidence is required for establishing the 
quantum nature of the object under study. Unlike this 
situation, with the superposition of three classically sepa- 
rated states in the superconducting ring, the characteris- 
tics of a radio-frequency superconducting qutrit detector 
(rf SQUTRID) will possess essential distinctions^, allow- 
ing one to state definitely their quantum origin. Par- 
ticularly, the quantum inductance extrema in the qutrit 
should be observed at external magnetic flux $e = n^o- 
Here we present experimental evidence and a theoreti- 
cal analysis for the fact that a rf SQUID with atomic-size 
ScS contact can be put into superposition of three dis- 
tinct states with a macroscopically large time of energy 
relaxation to lower levels and thus be turned into a rf 
SQUTRID. We explore the voltage-current and voltage- 
flux (signal) characteristics of this new device for the first 
time. Note that in the qutrit we study, the superposition 
states are formed due to the fast tunneling of fiux through 
potential barriers in the triple- well symmetrical potential 
in the phase space and the removal of the degeneracy of 
states of equal energy in each well, unlike in Ref.^, where 
the transmon with its three lowest energy levels was used 
as a qutrit whose superposition states were prepared by 
exciting the transmon from the base energy level to two 
higher ones. 

A substantial difference between the Josephson proper- 
ties of ScS and SIS contacts at low temperatures has been 
predictediSrJ^ in the microscopic theory of superconduct- 
ing weak links. The singular potential corresponds to 
the case of a "clean" ScS contact with a dimension d 
much smaller than the superconducting coherence length 
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Figure 1: (Color online) Circuit diagram of the SQUTRID. 
Experimental set-up is characterized by the following pa- 
rameters: Lt = 1.2 /iH, Ct ~ 630pF, Ip{t) — locostot, 
LOT /2tv = 5.794 MHz, Q = 302, M = k^/L^ = 1.52 nH, 
and Lq — 0.3 nH; the qutrit values Ic and C are discussed 
in the text. The inset shows the fully niobium 3D toroidal 
construction of the qutrit cell, where 1 is the pusher, 2 is 
the membrane, 3 is a part of the Lt coil, 4 is the toroidal 
quantization loop (cavity), 5 is the body, and 6 is the needle. 



^0 and the electron elastic mean free path £ (the so-called 
ballistic regime d <C ^). This fact leads to some pecu- 
liarities in the macroscopic quantum tunnelingi^ii^ and 
the quantum coherence of magnetic flux states in a su- 
perconducting ring closed by ScS contact. The major 
ones are an appreciable increase in the qutrit key per- 
formance parameter, the splitting of degenerate energy 
levels in separate potential wells and the emergence of 
high nonlinearity in the resulting qutrit superpositional 
levels that can be used in the rf SQUTRID based on rf 
SQUID circuitry (Fig. 1). 



II. THE SQUTRID MODEL 

Currently, two types of point contacts are distin- 
guished, depending on the ratio between the contact di- 
mension d and the electron wave length Xp = h/pp'- 
d ^ Xp for a classical point contact and d ^ Xp for 
a quantum point contact, In metals, actually, a quan- 
tum point contact is necessarily of atomic dimensions, as 
the electron wave length is of the same order of magni- 
tude as the atomic separation. For both classical^^ and 
quantumi^ ScS point contact with the critical current 1^, 
at r = the current-phase relation reads 



Is {(fi) = Ic sin ^ sgn[cos ^] , 4 



ttAq 
cRn ' 



where Aq is the superconducting energy gap and Rn is 
the normal-state resistance of ScS contact. The crit- 
ical current of the atomic-size (quantum) ScS contact 
was predicted^^^ to be quantized (as a consequence of 
the quantization of the contact conductance i?^^, in 
units of Go = 2e'^/h), Ic = iV(eAo//i), which was ob- 
served experimentally. From ^ one gets the Joseph- 
son coupling energy of ScS contact in the form Uj — 

-(/e$oA)|cOs((^/2)|. 

To develop the SQUTRID model in a zero-temperature 
approximation, we take the quantum Hamiltonian in the 
flux representation of the superconducting loop (Fig. 1) of 
inductance Lq closed by a clean atomic-size ScS contact 
with critical current Ic and self-capacitance C (so that 
parameter g = Ej/Ec = ^qIcC / {2ne'^) > 1) in the 
form^^-i^ 
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(2) 

where / = $/$o and fc — ^e/'&o are the normalized in- 
ternal magnetic flux $ in the loop and external magnetic 
flux <I>e applied to the loop. The quantum dynamical ob- 
servable of the internal magnetic flux in the loop is given 
by an operator of flux conjugated to an operator of charge 
in the contact capacitance: [$, Q] = —ih^ The key fea- 
ture of Hamiltonian ^ is its singular potential U{f\ fc) 
following from the nonsine current-phase relation ([T|) for 
ScS contact. Note that a model with both the potential 
attributed to ScS contact and the dissipation vanishing 
at zero temperature can satisfactorily describe the exper- 
iments on macroscopic quantum tunneling in a ring with 
a clean ScS contact, as shown in Ref. [ij. 

The solutions of the stationary Schrodinger equation 



Hq{f-fc)^{f)=E{fc)^{f) 



(3) 



(1) 



with Hamiltonian ^ yield wave functions 5'(/) and ener- 
gies E{fc) of the stationary states of the superconducting 
loop with ScS contact at a specified external magnetic 
flux fe- Let us consider the SQUTRID superconduct- 
ing loop where a three-well potential is formed. We get 
a series of states appeared during fast (with rf genera- 
tor rate to) increasing of external flux $e from to $0 
(Fig. 2) using parameters in Eq. ([3]) close to our exper- 
imental values: Lq = 0.3 nH, C = 4.36 fF, and = 4.0 
{Ic ~ AAfiA). It is seen from these solutions that, with 
a change in external flux, the initial state of a three- 
well symmetrical potential localized in the central well 
at $6 = [Fig. 2(a)] transforms through intermediate 
states [Fig. 2(b) and 2(c)] into a superposition state in a 
three-well symmetrical potential at $e = ^0 [Fig. 2(d)]. 
Note that the energy exchange rate between the two clas- 
sically separated states in a two- well symmetrical poten- 
tial at $e = 'i'o/2 is exponentially small at the specified 
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Figure 2: (Color online) Forming the base superposition 
SQUTRID state by rapid change of external flux "I>e- The 
set of nonstationary potential energies U{f)/kB, in terms of 
temperature, and squared wave functions are shown 

vs normalized internal flux / — $/$o at various $e; /9l = 4.0. 
Subscript after E refers to the energy level number, (a) The 
SQUTRID is initially in its ground state; its wave function is 
localized in the absolute minimum of the symmetrical poten- 
tial, (b) The potential becomes asymmetrical with two wells, 
but the system is still in its ground state, and the wave func- 
tion is localized in the minimum of the deeper well, (c) The 
wave function remains surprisingly localized in the same sin- 
gle well since the potential changes too quickly in comparison 
with the tunneling rate (the interwell barrier is rather high in 
this case), (d) The base state |^'so|^ of the three-well super- 
position is formed in the potential symmetry point when the 
barriers are small, and tunneling time is consequently short. 
The relaxation time of this state turns out to be macroscop- 
ically large, so that it becomes stable enough and "latched" 
for many further cycles of <E>e. 



parameters, so the system state remains localized in the 
starting potential well during the increase in external flux 
toward the point $e = $o- In this point, the superpo- 
sition qutrit state |*so) = ci|*i) -l-C2|*2> + cal^s) of 
the wave functions of all three separate wells is formed in 
the three-well symmetrical potential, similar to the for- 
mation of the superposition qubit state in the two-well 
potential. 

Quantum coherence of the qutrit flux state |4'so) in 
a three-well potential manifests itself as coherent oscil- 
lations of magnetic flux between all the three potential 
wells due to its fast tunneling through the potential bar- 
riers separating the central and the side wells. Numerical 
analysis of Eq. ([3]) shows that, for resonant tunneling in 
a three-well potential of the superconducting loop of in- 
ductance Lg = 0.3 closed by a clean ScS contact with 
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Figure 3: (Color online) Superposition of three states in the 
superconducting loop closed by a clean ScS contact. Distri- 
butions of squared wave functions |*I'so(/)|^ among the three 
potential wells vs normalized internal flux / = $/$o for two 
close values of external magnetic fluxes $e corresponding to 
(a) fully symmetrical and (b) slightly tilted potential U{f); 
Pl = 4.0. Superposition energy levels are denoted Eso,Esi, 
and Es2- The process of energy relaxation of the base su- 
perposition state ^so(/) from level Eso to level E-i with 
characteristic time is shown by arrows, (c) Superposition 
energy levels Eso,Esi, and Es2 vs normalized external mag- 
netic flux fe = "I'e/'&o- (d) Normalized quantum inductance 
of the base superposition level Eso vs fe : curve 1 is zero- noise 
(LqLq^){fe) and curve 2 is noise-affected {LqLq^)ef f{fe) av- 
eraged over low-frequency noise with a = \/{5fi) = 0.013. 



C = 3 — 6 fF {/3l = 4.0), the flux oscillation rate between 
wells I'Qi = AEoi/h reaches 25 — 40 GHz. The magnetic 
moment of the flux states corresponding to the side wells 
of the three-well potential 4^'^^ = I^S ~ 10"" J/T~ 
10^^ fj,B (where Is ~ 0.6$o/ig ~ 4/LtA is the supercur- 
rent in the side-well flux states, S* ~ 2 x 10^^ m^ is the 
loop area, and fiB = 0.93 x 10^^'^ J/T is the Bohr mag- 
neton), and the magnetic moment of the central flux 

state /is = (since the supercurrent is zero in this 
state). Thus, we have the coherent superposition of wave 
functions corresponding to the three distinct macroscopic 
flux states in the three- well symmetrical potential of the 
qutrit at $e — n^Q. 

Let us refer to |^so), which is the three- well su- 
perposition state with minimum energy, as the "base" 
SQUTRID superposition state since adiabatic movement 
along the respective energy level Esq plays the main role 
in our experiments. State l^'so) in the potential U{f; /e) 
is characterized by quantum number n ^ 1 [n = 30 
for the above-cited parameter values; see Fig. 2(d)]; i.e.. 
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the deep central well {AU/kB — 300 K) contains a large 
number of quantum levels. The time of the energy re- 
laxation of the base SQUTRID state to underlying states 
in the central well must be macroscopically large to en- 
able measuring the superpositional nonlinearity. This is 
achieved due to the design of the qutrit loop in the form 
of a high-quality three-dimensional (3D) toroidal super- 
conducting cavity (see inset in Fig. 1 and Sec. Ill), which 
has no resonant modes with frequencies corresponding 
to the frequencies (energies) of transitions from the base 
superpositional level Esq to underlying energy levels. 

Figure 3(a) displays numeric solutions of Eq. ^ for the 
squared absolute value of the wave function of the base 
superposition state |4'so(/)P and energy levels £^50: Esi, 
and Es2 of the three superposition states as well. The 
fact that splitting AEoi/ks = [Esi - Eso)/kB = 1.65K 
between the base and the nearest superposition energy 
levels substantially exceeds the environment tempera- 
ture for the chosen parameters is key in the present 
experiment {AEoi/kBT = 1.65/0.45 « 3.7). As a re- 
sult, the broadening of level Esq because of statisti- 
cal averaging over equilibrium density matrix ^£'50 — 
A£;oi/[l + exp{AEni/kBT)] < AEqi is relatively smah 
{SEso/kB ~ 0.04 K) and can be neglected in the analy- 
sis of the system energy structure (zero-temperature ap- 
proximation). Also, it is obvious that the large tunnel 
splitting, which multiply overruns the temperature value, 
guarantees good "quantumness" of the superconducting 
loop closed by a clean ScS contact provided that the en- 
vironment and the measurement circuit noises are effec- 
tively suppressed. At the same time, the base superposi- 
tion level Esq is situated far below the potential barrier 
top Ub, {Ub — Eso)/kB ~ 5K T, so that thermal tran- 
sitions rate over the barrier is vanishingly small, and the 
process of relaxation of the metastable base superposition 
state into the deep central well can be strongly limited 
for the chosen parameters since {Esq — E^ij/ks « 8K. 

Small adiabatic variations of the external flux relative 
to symmetry points $e = lead to tilting the three- 
well potential, with a change in the wave functions and 
the energy levels. As seen from Fig. 3b, when the flux 
shifts apart from the symmetry point $e = 'I'oi the wave 
packet of the base state |^so(/)P partially transfers into 
the deeper side potential well. A rearrangement of the 
superposition wave packets and thus the formation of the 
superposition energy levels Esoife), Esi{fe) and £'s2(/e) 
[shown in Fig. 3(c)] occur over the time span reciprocal 
to the flux tunneling rate, ts ~ i'qi ■ In experiment, we 
deal with relatively low-frequency (adiabatic) processes, 
in which the time-dependent superposition levels practi- 
cally coincide with the stationary ones. 

The superconducting qutrit detector is built on the 
basis of rf SQUID circuitry. That is, the qutrit is in- 
ductively coupled to a high-quality LtCt superconduct- 
ing tank (Fig. 1) serving as a linear classical detector 
{hojT <C ksT). Only a small fraction of the flux in the 
LtCt tank (M/Lt ~ 0.0013 in our experiment, where 
Lt is the tank inductance and M is the mutual induc- 



tance between the qutrit and the tank) is transferred to 
the qutrit, which allows making use of the concept of 
weak continuous quantum measurements i^^rJ^ When the 
characteristic frequency of the classical detector {LtCt 
tank) is much lower than that of the quantum object 
(in our system the ratio uj /vqi ^ 10^^), the latter's dy- 
namics can be studied by means of quantum-mechanical 
equations with the detector classical field treated as the 
external adiabatic parameter)^"— Within this semiclas- 
sical approach of treating the qutrit plus tank system, 
one obtains the classical equation for the LtCt tank 
containing the parametric quantum inductance contribu- 
tion [see Eq. (O], instead of the parametric Josephson in- 
ductance contribution probed in rf SQUID by means of 
the impedance measurement technique (IMT). The IMT 
with a weak continuous quantum readout was success- 
fully applied to studying different types of qubits.— 

For the LtCt tank driven by rf-bias current Iq cos cut, 
the measured output functions are the amplitude of volt- 
age oscillations Vt and the voltage-current phase shift 
aT- The equation for voltage V{t) across the LtCt tank 
having the quality factor Q = ujtRtCt, with the current 
contribution due to the weakly coupled qutrit Mlqo/LT 
taken into account, reads 

CtV+-^ + -^ I Vdt+^^ - Ip{t), Ip[t) = locosLot. 
Kt Lt J Lt 

(4) 

Here 

r d{^so\H,\^so) _ dEso{'^>e) ... 

^.0(<fe) - - (5) 

is the current circulating in the superconducting loop of 
the qutrit in its base quantum superposition state, as a 
function of the total external flux $e. 

In further considerations it is convenient to use the 
function ig^(/e), which is called the reciprocal quantum 
inductance, defined as 

(6) 

This function, being, in fact, the local curvature of energy 
level Esoife), describes the nonlinear properties of the 
qutrit in the base quantum superposition state at small 
variations of the external magnetic flux. At the same 
time, it characterizes the parametric inductance inserted 
in the LtCt tank due to a weakly coupled quantum de- 
vice. 

Considering the emf induced in the qutrit loop $e = 
MIl = MV{t)/LT, we obtain 

v+uj^tV = fiy,v,t), 



f(y, V, t) = -e L,LQ\'f,{t)] Lo^V - 




M f 

$e(t) - + $ac(i) = + — / Vit)dt , 

(7) 
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where is the direct (low-frequency signal) external 
flux biasing the qutrit loop and $ac(^) is the alternat- 
ing external flux applied to the loop due to the tank 
flux oscillations. Thus, the strongly nonlinear reciprocal 
quantum inductance function LQ^(/e) [Eq. ([6])], charac- 
terizing the curvature of the qutrit base superposition 
energy level, will determine the solution of Eq. ([T]) . If 
the condition f{V,V,t) ^ uj'^V is fulfilled, which is 
valid when Q > 1, /q < uitCtVt, k'^LqLg^ < 1, the 
Krylov-Bogolubov method for solving weakly nonlinear 
equations^^ can be applied to solve Eq. (O . Substituting 
voltage V{t) in the form 



dynamic fluctuations.^^ In the Gauss-distributed noise 
approximation we get 



V{t) = Vrit) cos[wt + arit)], 



(8) 



where Vr (t) and (t) are slowly varying functions (with 
small relative variation over the oscillation period T — 
2tt/ujt), we get abridged equations for Vr(i), Oirit), and 
the equations for the voltage amplitude and phase shift 
of steady-state [Vrit] — 0, axit) — 0] oscillations in the 
LtCt tank: 



Vt 



(jjtLtQ Iq 



v/l+4e(yT,$<ic)2Q2 



2ir 




dc 



tan ax 



MVt 
ojLt 



cos^ T dr. 



(9) 

where Co = (1 ~ /^) is a detuning parameter set to zero 
hereinafter since uj ~ ujt in experiment. As seen from 
Eqs. dH), voltage-current Vt{Io) and voltage- flux (signal) 
Vri^dc) characteristics of the SQUTRID are determined 
by the reciprocal quantum inductance LQ^($e) averaged 
over a period of oscillations in the tank. Due to the sharp 
dependence of L'^i^e) in the vicinity of = [see 
Fig. 3(d)], small variations of signal magnetic flux 
will lead to a substantial change in the reactive part of 
the tank impedance and therefore in the Vt (t) and ax (t) 
dependencies. 

Equations ([9|) should be solved numerically because of 
the strong nonlinearity of the LQ^{^e) function, to which 
the sought tank voltage amplitude Vt enters through 
the external flux <I>e. We are also interested in taking 
into account the effect of noise (generally of complex na- 
ture and spectrum) influencing the qutrit on the mea- 
sured Vt{Io) and Vri^dc) dependencies. To this end, 
a simplified model is used in which the major part of 
the noise influencing the qutrit loop is considered to be 
caused by the measurement circuit. The noise from the 
circuit produces fluctuations of external flux applied to 
the qutrit loop that change the qutrit quantum response 
and, in turn, its back action to the LtCt tank. If 
the inequality lut ^ ^ AEqi/H is valid for all the 
noise spectrum components oJi affecting the SQUTRID 
loop from the side of the measurement circuit, then the 
task becomes easier, and the value (LgLQ^)e//(/e) ef- 
fectively contributing to the tank can be found using 
the method of averaging over quasistationary thermo- 



(^9-^Q^)e//(/e) 



(10) 

where a = y/ {Sf^) is the standard deviation of the 
noise flux associated with the measurement circuit (and 
is the noise flux variance). Substituting this func- 
tion, {LqLQ^)ef fife), of the effective reciprocal normed 
to Lq quantum inductance [see Fig. 3(d)] into Eq. ([S]), 
one can numerically obtain voltage-current Vt{Io) and 
voltage- flux Vri^dc) characteristics of the SQUTRID 
that account for the measurement-induced noise flux 
parametrized by its variance. 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 



The main objective of the experimental design is to 
provide conditions at which (i) degenerate levels of the 
separated wells become split with AEqi ksT due to 
interwell tunneling and (ii) all the frequencies in the spec- 
trum of the environment noise Ui will be small compared 
to the rate of transitions to the superposition level Esi 
and the lower level E^i situated in the middle well, i.e., 
uji <C AEoi/h; {Esq — E-i)/h. In order to meet these 
conditions, a 3D toroidal SQUTRID loopi^ with induc- 
tance Lq — 0.3 nH was made from pure (99.999%) nio- 
bium (Nb) . The design of the 3D toroidal superconduct- 
ing loop (inset in Fig. 1), which, in fact, becomes a closed 
3D cavity, favorably eliminates undesirable coupling of 
the qutrit to the external electromagnetic environment. 
The only way for electromagnetic interference to come in 
is through the narrow (0.5 mm) and long (8 mm) chan- 
nel made in the qutrit body for wiring the coupling coil 
placed inside the qutrit toroidal cavity. With these di- 
mensions, such a channel acts as a below-cutoff waveg- 
uide for all frequencies of up to several hundred gigahertz 
providing very high attenuation. The coupling coil leads 
are properly filtered. Adjusted in situ atomic-size ScS 
contact was realized between a surface-cleaned Nb nee- 
dle and an annealed Nb plate with a crystallite size close 
to 0.5 mm. A number of atoms N in the contact opening 
can be assessed as ratio of the contact critical current 
Ic to its quantizing value (eAo/fi), giving TV ~ 10. The 
sample, the main part of the resonance tank, and the fil- 
ter were cooled down to T = 450 mK in a '^He pumped 
refrigerator cell (see Fig. 1). The first stage of the rf am- 
plifier and additional filters were placed at T = 1.5 K. A 
triple /i-metal shield around the liquid "^He Dewar and 
a superconducting Pb shield around the measuring cell 
were used to reduce and stabilize the ambient magnetic 
field. 

An unusually low resonance frequency w^/Ztt — 
5.79 MHz of the LtCt tank was chosen to decrease the 
potential variation rate and to meet the adiabatic con- 
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editions for the qutrit quantum dynamics even at high 
(such that MIqQ ~ $o) amphtudes of the rf generator 
current Iq. The calculated time (t ^ 10~^s) of passing 
the dip of the function {LqLQ^)ef f{fe) [Fig. 3(d)] for this 
frequency considerably exceeds the superposition setting 
time (rs ^ 3x 10^^^ s) in the three-well symmetrical po- 
tential [Fig. 3(a)]. The resonance tank permanently mea- 
sures the state of the quantum system which is weakly 
coupled to it and, at the same time, generates additional 
noise in the superconducting loop. Thus, we should ex- 
pect that the measurement process limits the quantum 
superposition in our system. However, detecting the av- 
eraged curvature of the base energy level is still possible 
since the uncertainty of the magnetic flux associated with 
the effect of the measurement circuit and temperature is 
estimated to be as low as ^10^^$o. 

Figure. 4(a) exhibits a set of experimental voltage- flux 
characteristics Vri^dc) of the SQUTRID with /?l ~ 4 
obtained while sweeping the external magnetic flux 
for several rf generator current amplitudes Iq . Note that 
plateaus exist in the range around ^dc = 'i'o("- + 1/2), 
shrinking with increasing Iq- These plateaus corre- 
spond to a quasiautonomous LtCt tank with charac- 
teristic resonance impedance Rc — Vt/Iq — 13.3 /cfi 
and indicate that no inductance is inserted into the 
tank from the qutrit that receive external flux $e(i) = 
<^dc + <^ac{t), $ac(i) = MV{t)/{uJTLT), while Sweep- 
ing ^dc flux. This, in turn, is evidence of the function 
{LqLQ^)eff{fe) becoming zero [see Fig. 3(d), curve2] and 
of a no-superposition qutrit state in the respective range 
of ^e{t). The onset of nonlinearity in the Vri^dc) curve 
is observed at ^e{t) = ^dc+^ac{t) ~ n^o, when the mag- 
netic flux <&e(i) falls into the narrow region around 7i$o 
[Fig. 3(d), curve 2] where the qutrit superposition nonlin- 
earity is localized. Particularly, for the triangle-shaped 
signal curve [at Iq = 2.55 nA in Fig. 4(a)] with a plateau 
degenerated into the point {^dc = 0.5<f>o, = 30/xV), we 
have $e = ^dc+^ac « ^-o, where <^ac = MVt/{ujtLt) ~ 
0.5$o is the amplitude of the ac flux coming to the qutrit 
from the tank. As clearly seen in Fig. 4a, the Vri^dc) 
characteristics of the tank are maximally affected by the 
quantum inductance inserted from the qutrit at symme- 
try points ^dc = n^O: where the qutrit superposition 
nonlinearity is maximum. At low-to-moderate current 
amplitudes, Vri^dc) dependencies are well described by 
the theoretical model [see Fig. 4(b)] that takes into ac- 
count the noise influence of the measuring channel, with 
independently measured SQUTRID parameters and ca- 
pacitance C = 4.36 fF being used. 

Figure 4(c) presents the initial parts of the SQUTRID 
Vt{Io) characteristics (/3l w 4) registered for two values 
of the magnetic flux, <^dc = 'i'o and ^dc = ^of^- The 
theoretical curves Vt{Io) are derived from Eqs. ([9|) with 
experimentally measured SQUTRID parameters and av- 
eraging the superposition nonlinearity of the quantum 
system over low-frequency noise ()10|) with standard de- 
viations a = 0.013 and a = 0.023. For the V-I curves 
with ^dc = n^o, the effective reciprocal quantum indue- 
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Figure 4: (Color online) SQUTRID characteristics: (a) exper- 
imental voltage-flux (signal) characteristics Vri^dc/^o) ob- 
tained for several amplitudes of rf generator current Jo (indi- 
cated as the curve parameter) and (b) set of voltage-flux char- 
acteristics calculated accounting for noise with standard devi- 
ation a — 0.013 at currents Iq close to the experimental values 
(curve parameter). Independently measured SQUTRID pa- 
rameters Lg = 3x 10"^° H, Pl = 4.0, Lt = 1.2 x lO"** H, 
M = 1.52 X 10"^ H, g = 302, fc^Q = 1.89, and capacitance 
C = 4.36 fP were used in the calculations, (c) Experimental 
(dots) and theoretical voltage-current characteristics Vt{Io) 
calculated for noise flux standard deviations a = 0.013 (dot- 
ted line) and a = 0.023 (dashed line) at $dc ~ $o and no- 
superposition experimental (squares) voltage-current charac- 
teristic Vt{Io) at $dc ~ "l'o/2; (d) smearing of Vt{Io) curves 
when the refrigerator temperature rises to T = 1.05 K. Non- 
linearities due to the superposition (circled) almost vanish. 



tance of the base superposition qutrit level, introduced 
into the tank at low rf generator currents {Iq < 0.4 nA), 
leads to a large shift in the tank resonance frequency 
from the generator frequency, which results in a de- 
crease in the voltage detected at the LtCt tank down 
to a value comparable to the noise level. As seen from 
Eqs. (HI) and pI7|) . at ^dc = n<^o and a low cur- 
rent (linear) regime the voltage reads approximately as 
Vt « LJTLTloQ/lk^QiLqLQ^ffl; that is, it is reduced 
by a factor of \PQ{LqLQ^)^^^\ > 1 compared to the 
autonomous LtCt tank voltage. 

The theoretical curve Vt{Iq) with noise standard de- 
viation cr — 0.013 coincides with the experimental curve 
in the range of low-to-moderate generator currents. At 
higher a, the factor \k'^Q{LqLQ^)^j^\ decreases, and the 
slope of the V-I curve in the area of small rf genera- 
tor currents becomes higher than in experiment [see the 
curve with a = 0.023 in Fig. 4(c)]. At Iq > 2.5nA, the 
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amplitude of ac magnetic flux induced in the SQUTRID 
loop $ac — MVt[Io) / {o-itLt) exceeds the half-width of 
the "main dip" of the effective reciprocal quantum in- 
ductance, $oc > 0.02$o [see Fig. 3(d)], and then Vt 
rapidly rises with the increase in the rf generator current. 
Such behavior at high rf generator currents /q (nonlin- 
ear regime) is well described by the proposed theoreti- 
cal model. The observed quantitative deviation can be 
attributed to both simplifications made when deriving 
the theoretical model (ideal ScS contact, zero SQUTRID 
temperature, and Gaussian noise distribution) and an ex- 
perimental error in determining the SQUTRID parame- 
ters. The no-superposition quasiautonomous-tank exper- 
imental V-I curve at <^dc = $o/2, which is related to the 
considered "plateau regime" in Vri^dc) characteristics, 
is shown in Fig. 4(c). It is seen that this V-I curve follows 
Ohm's law at the specified /3l ~ 4. As discussed above, 
with a high enough barrier, separating the two wells in 
the two-well potential, the qutrit wave function remains 
localized in one of them during the variation of the po- 
tential by the ac component of the external magnetic 
flux. 

The impact of thermodynamic fluctuations on the ef- 
fective reciprocal quantum inductance (|10p and on the 
rf V-I characteristics of the SQUTRID is illustrated in 
Fig. 4(d). The nonlinearity due to the superposition of 
states is severely smeared out by temperature at T = 
1.05 K and becomes completely unobservable in our ex- 
periments at helium bath temperature T > 1.5 K. This 
result can be explained by an increase in the higher- 
level population because of the reduction of the ratio 
AEoi/kBT < 1 and an essential increase in the noise 
flux variance a^. 

To see nonlinearities caused by the superposition of 
states in more complicated potentials, we have increased 
the critical current of the SQUTRID ScS contact up 
to Ic ~ 8 /iA. An increased Ic and, thus param- 
eter corresponds to contacts with a greater number of 
atoms (atomic rows) at their opening and hence de- 
creased values of their normal-state resistance Rn- Fig- 
ure 5(a) displays rf V-I characteristics of the SQUTRID 
with /3i w 7.3 obtained at temperature T = 450 mK for 
magnetic fluxes $dc = n^o and $dc = {n +1/2) $0- The 
initial parts of both branches of the V-I curves remain al- 
most linear due to dynamic localization of the qutrit wave 
function. The effect of the reciprocal quantum induc- 
tance determined by the effective curvature ([TO)) of the 
base superposition level Eso is observable at generator 
currents less than those needed to form "classical" steps 
in SQUID rf V-I characteristics. Figure 5(b) shows super- 
position energy levels calculated for /Sl = 7.3 that appear 
in the four- well symmetrical potential at ^e{t) — 1.5$o- 
Despite the fact that, at chosen SQUTRID parameters, 
the environment needs to absorb a considerable amount 
of energy [{Eso — E-{)/kB ~ 7K] when the system re- 
laxes to the lower level, the relaxation events often oc- 
curred in the experiment. This can be caused by very 
high numbers of superposition energy levels (n = 100 in 




Figure 5: (Color online) Superposition of states in the su- 
perconducting ring closed by ScS contact, with a four- well 
potential, (a) Experimental voltage-current characteristics of 
the SQUTRID with ~ 7.3 obtained at T = 450 mK for ex- 
ternal magnetic fluxes $dc = $o/2 and ^dc ~ '&o- The areas 
where the eff'ect associated with maximal local curvature of 
the base superposition level appears are circled. The jumps 
shown by the arrows are due to relaxation events, (b) Poten- 
tial energy and superposition energy levels calculated for the 
SQUTRID parameters taken from experimental data. Three 
superposition levels still form in this four-well potential and 
are situated far above two middle wells. 



our case) that are less stable with respect to noise and by 
possible dissipation effects related to an increase in the 
number of atoms in the contact. The effect of the energy 
relaxation resulting in a jump to a linear part of another 
{^dc ~ 0) no-superposition branch of the V-I curve is 
shown in Fig. 5(a) by the up arrow. The reverse process 
of the transition from the no-superposition V-I branch to 
the superposition one shown by the down arrow reflects 
the formation of the superposition state in the four-well 
symmetrical potential conflguration at a driving current 
that gives the total external flux ^e{t) — 1.5$o- 

Note that small amplitudes of the SQUTRID driving 
current in quantum measurements are the most favorable 
for quantum informatics. The obtained conversion ratio 
{j^ = dVT/d^e « 3x10^° V/Wb) and the SQUTRID per- 
formance in the three-well superposition (at (3^ ~ 4) can 
be further increased by increasing the frequency uj while 
simultaneously and proportionally reducing the gener- 
ator current /q. For example, at pumping frequency 
w/27r = 100 MHz we have 77 of up to 5 x 10"V/Wb. 
In this case the SQUTRID voltage-flux characteristics 
would be similar to those shown in Figs. 4(a) and 4(b) 
for small (Iq < 0.4 nA) current amplitudes. Since, as 
shown in Ref.[2^ the detector and the tank temperature 
in these measurements may be as low as 10 mK and the 
SQUTRID acts as a sensor of parametric quantum in- 
ductance under adiabatic conditions, such a device can be 
considered to be a representative of the class of quantum- 
limited detectors. 
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IV. SUMMARY AND CONCLUSIONS 

We have constructed and studied a new ideal paramet- 
ric detector of magnetic flux named rf SQUTRID which 
is based on quantum superposition of three macroscopic 
flux states of a 3D toroidal superconducting loop closed 
by Nb-Nb atomic-size point contact. Due to the specific 
form of the potential barrier, the small contact capaci- 
tance, and the fast tunneling dynamics of formation of a 
coherent three- well superposition state {ts ^ 3xl0~^^s), 
the adiabatic conditions in the rf SQUTRID are valid 
up to frequencies u}/2tt « 3 GHz for small driving am- 
plitudes (M/o ^ 10"^$o)- Therefore, the flux detec- 
tors based on the dependence of the local curvature of 
the base superposition energy level (or quantum induc- 
tance) on the external magnetic flux at low temperatures 
fall well within the class of fast and sensitive devices. 
Note that unlike the rf SQUTRID/SQUBID based on the 
atomic-size point contact, similar detectors with contacts 
of the SIS type will have slow tunneling dynamics of flux 
wave packets and, respectively, a lower operating speed. 
The rf SQUTRID (as well as rf SQUBID) is the flux 
magnetometer "dual" to the rf single electron transistor 
(SET) electrometer;^ and so the main applications of 
the rf SQUTRID will initially be experimental study of 
physical processes in flux qubits and engineering readout 
systems (sensors) for reading states in the small-scale su- 
perconducting quantum registers using weak continuous 
measurements. 



We stress the following results of this paper: (i) The co- 
herent superposition state 4*50 (/) of wave functions cor- 
responding to three distinct macroscopic states is formed 
in the three-well symmetrical potential conflguration at 
$e ~ in the superconducting circuit with the clean 
atomic-size ScS contact under study, (ii) The energy re- 
laxation time Te of the base superposition state ^'so(/) 
to lower energy states can be macroscopically large to 
perform quantum measurements, (iii) Unlike classical rf 
SQUID in a nonhysteretic regime, the nonlinear proper- 
ties of the superposition in the three-well potential un- 
der the condition cj/27r <C vqi ^ Aq/Zi enable making 
a sensitive and fast parametric detector of magnetic flux 
without a quasiparticle current, i.e., principally a nondis- 
sipativc rf SQUTRID device, (iv) An unusually large 
splitting AEoi/kB — 1.5 — 2K and a very large (at low 
noise variance) LqLg^ « 200, specifying the nonlinearity 
of the studied quantum system, make the rf SQUTRID 
with a ScS contact a very promising element for quantum 
informatics. 
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